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The  ever  increasing  need  to  be  informed  has  caused  an 
ever  increasing  volume  of  publications.  The  very  quantity 
of  these  publications,  however,  makes  it  difficult  to  satis¬ 
fy  the  need,  because  it  becomes  more  and  more  arduous 
to  locate  the  literature  pertinent  to  any  one  subject;  and 
because  any  one  of  us  becomes  more  and  more  weary  of 
having  so  much  to  read.  Moreover,  only  a  small  minority 
of  these  publications  may  help  satisfy  the  need,  because 
an  overwhelming  majority  is  void  of  content. 

Consequently,  it  is  imperative  to  present  to  the  public 
only  manuscripts  which  do  carry  new  information  without 
undue  verbiage;  hoping  to  be  judged  not  by  the  number 
of  published  pages  but  by  their  quality.  In  particular, 
papers  in  the  field  of  Applied  Mathematics  ought  to  be 
published  only  if  they  contain  one  or  more  of  the  follow¬ 
ing  items:  new  basic  results,  new  methods,  new  applicat¬ 
ions,  new  numerical  results,  new  presentation  of  difficult 
and  important  topics,  up-to-date  bibliographies;  and  if 
the  number  of  their  pages  is  not  dictated  by  the  desire  of 
imposing  upon  the  superficial  reader. 

To  discharge  our  contractual  obligations,  we  publish 
Technical  or  Scientific  Reports,  such  as  the  one  you  now 
have  in  your  hands.  It  has  been  our  constant  policy  to 
see  to  it,  that  they  satisfy  the  above  strict  criterion. 
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Mathematical  Methods  in  Transducer  Field  Theory: 
The  Finite  Cylinder 

by 

N.  G.  Parke  III  and  W.  Williams 


Introduction 

This  report  deals  primarily  with  a  method  of  determining  the  velocity 
potential  for  the  sound  field  generated  by  a  finite  circular  cylinder  with 
rigid  end  faces  which  is  executing  pure  radial,  harmonic  vibrations  uni¬ 
formly  over  its  cylindrical  surface;  it  will  be  concluded,  following  the 
presentation  of  a  specific  numerical  calculation,  with  a  simple  extension 
of  the  method  applicable  to  the  determination  of  the  velocity  potential 
for  the  sound  field  generated  by  arbitrary  finite  bodies  of  revolution  vibrat¬ 
ing  harmonically  in  time. 


Background 

As  a  consequence  of  specifying  harmonic  time  dependence,  the  problem 
of  determining  the  time  dependent  velocity  potential  Y}i)  in  a  domain 
exterior  to  the  vibrating  source  from  the  wave  equation: 

(1)  -  o 

(in  which  e  is  the  sound  speed)  reduces,  upon  the  usual  substitution  of 

(2)  £  at)  -  [*<*;<? -‘■•*7 


in  (l)  to  the  determination  of  the  time  independent  velocity  potential  y>(?) 
in  a  domain  exterior  to  the  vibrating  source  from  the  Helmholtz  equation: 


(3)  +  HcX  p(r)  -  O 

(in  which  the  wave  number  $  =  where  A  is  "the”  wave  length  of  the 

sonic  radiation).  Further,  if  the  Sommerfeld  radiation  condition 


(4) 


r  w(t) 

\im  r  L  -3 f- 


o 

> 
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is  Imposed  upon  y>(r)  then  an  important  paper  by  Hartman  and  Wilcox' 
shows  that  tor  a  space  of  three  dimensions  the  following  series  expansion 
converges  in  the  mean  to  yiC?)  in  an  exterior  domain 

O-*  )n  « 

(5)  ^(r)  ~  2  2  iJ  !nC>(tr)  *  y  (tt  *) 

AC** 

(in  which  h.  is  a  spherical  Hankel  function  of  the  first  kind  and  V  is 
a  spherical  harmonic). 

Hence,  the  problem  of  determining  the  time  independent  velocity  potential 
tp(r)  in  domain  exterior  to  a  finite  source  of  the  sound  field  consists  in 
the  specification  of  the  constants  Cl.  subject  to  the  application  of  appropri- 
ate  boundary  conditions. 

The  Finite  Cyl:  -^er  Problem 

General  Theoretical  Analysis  -  In  order  to  specify  the  precise  boundary 
conditions,  the  center  of  the  circular  cylinder  of  radius  &  and  length  J.l> 
will  be  placed  at  the  origin  of  a  coordinate  system  in  which  either  the  cylin¬ 
drical  coordinate  set  &t>cp  or  the  spherical  coordinate  set  0,  will 
be  employed  depending  on  efficacy.  With  the  assertion  of  rigid  (i.e.,  motion¬ 
less)  end  faces  and  pure  radial  harmonic  vibratory  motion  of  the  cylindrical 
surface,  the  boundary  conditions  on  the  finite  cylinder,  which  follow  directly 

PS 

from  the  relation  between  the  velocity  v  and  the  time  dependent  velocity 
potential  (  ?,  ~t )3  v  t)  >  can  be  written,  apart  from  a  scale  factor, 


as 

(6a) 

dV'(r)  =  o 

dt 

on  the  end  surfaces:  O  -  f>  —  d.  3  2=*^ 

H 

(6b) 

dtG)  „  ) 

on  the  cylindrical  surface:  f  *  d.  ■  -  *•  L 

for  the  time  independent  velocity  potential  ^(rJ. 


*'ph.  Hartman  and  C.  Wilcox,  "On  Solutions  of  the  Helmholtz  Equation  in 
Exterior  Domains",  Math.  Zeitschr.,  75 >  228-255  ( 1 9^1 ) . 
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Two  immediate  simplifications  in  the  structure  of  the  series  expansion 

(5)  for  follow  immediately  from  symmetry  considerations: 

a)  The  specification  of  the  cylinder  as  circular  dictates  azimuthal 
symmetry;  consequently,  wi«o  in  (5)  and 

where  ?/*)  is  the  usual  Legendre  polynomial. 

b)  The  form  of  the  Neumann  condition  (6b)  on  the  cylindrical  surface 
dictates  equatorial  symmetry  (i.e.,  symmetry  around  the  plane 

O -  t  a 7 r  or,  equivalently,  the  plane  t  O-  *tr  ); 

consequently,  (f,  -i  )  — »  $COu*)  -  ^  . 

Hence, 

(7)  [£>( r)  »  <p  (rtb)  At  CCts  dJ  eve* 

(in  which  (*■  has  been  absorbed  into  dh  and  where  hn(1*r)  without  the 
super  script  (l)  will  continue  to  designate  the  spherical  Hankel  function 
of  the  first  kind),  and  the  problem  for  the  finite  cylinder  is  reduced  to 
the  two  dimensional  problem  -  because  of  the  azimuthal  symmetry  -  of  finding 
the  complex  coefficients  2#)  «t4 ,  by.,...,  which  satisfy  the  Neumann  conditions 

(6)  on  the  rectangular,  cross  sectional  boundary  of  the  cylinder. 

In  order  to  facilitate  the  analytical  development  which  is  required  to 
calculate  the  normal  derivative  of  the  series  expansion  of  (7),  evaluate  it 
on  the  rectangular  boundary  and  match  this  evaluated  derivative  to  the  "data" 
for  the  problem  as  given  by  the  boundary  conditions  of  (6),  the  polar  equation 
of  the  rectangular  boundary  is  derived  as  follows.  Examination  of  Figure  1 
shows  that  the  rectangular  boundary  -  and  also  the  domain  exterior  to  the 
cylinder  -  is  naturally  divided  into  three  regions: 

Region  I:  O  —  d  —  ^ 

(8)  Region  II:  »  a  d  ^  77-1%, 

Region  III:  Tf-6,  Z  7T 
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Figure  1  -  Cross  Section  of  Finite  Cylinder 


(in  which  the  polar  angle  6  is  measured  from  the  positive  2  -axis)  and 
that  the  radial  distance  from  the  coordinate  origin  to  the  boundary  of 
the  rectangular  cross  section  is  given  in  theae  separate  regions  by: 


(9) 


Region  I: 

r  .  =  U/cox a 

«y/  > 

Region  II: 

rey|  «=  a  , 

Region  III: 

rcyi  -  "W 

where  the  negative  sign  is  introduced  in  Region  III  in  order  to  maintain  Y?o 

in  6s  second  quadrant.  If  a  discontinuous  function,  designated  by  us  as  an 

(2) 

H-function' , '  is  introduced  with  the  property  for  the  variable  X  that 


(10) 


h  ( r} ». ) 


Jfi  x  i  h 
oUerwitc 


then  the  three  polar  equations  for  the  rectangular  boundary  of  (9)  can  be 
combined  into  the  form 

(11)  l/coj4  Lh<#,0  ”  HC-rr-a.,  ir)]  +  a/s;»a  W  (t,  7 T- e*  ) 


W) 


Which  is  the  characteristic  function  of  a  set. 
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(in  which  0.  *  €  where  €  is  infinitesimally  small^^).  If  the 

notation 

H(o,  S.)  — H.Ca)  — )  H, 

H(«.f  TT-d)  —7  H4(aJ  —7 
HCtt- a#)  it)  Hs 

is  employed,  then  the  polar  equation  for  the  boundary  becomes 

( 12 )  f  y/  — ^  3  =  &  (  H,  -  )  +  */s;n  o 

which  is  the  desired  result. 

In  addition,  the  use  of  the  H-functions  leads  immediately  to  the 
following  expressions  for  p  and  i  on  the  rectangular  cross  section 

(13a)  j*  =  b  +*»6  (H,  -  Wj  )  -b  S  Hj, 

(13b)  2=  U  H.-rtJ  +  S  c»te  H4 

which,  together  with  the  obvious  inter-regional  property  for  the  H-functions 

(1M  Hi  Hi  -  sis  Hi 

(in  which  S-j  is  the  Kronecker  delta),  leads  to  the  relation 

(15)  S2-  b  / 2«  a  Ch,  -r  iV*  )  +  d.  /S/t>i&  -3, 

which  is,  of  course,  also  directly  obtainable  from  (12). 

The  normal  derivative  ,  which  will  be  calculated  over  the  boundary 

and  then  matched  to  the  boundary  conditions,  has  the  form 

06)  sj-M  =  (  +  K 

3  H  31 

(in  which  the  negative  sign  is  introduced  to  force  the  normal  in  the  negative 
i  direction  in  Region  III).  If,  now,  and  are  applied  directly 


The  distinction  between  &e  and  6m  has  been  emphasized  so  that  the  regional 
definitions  of  (8)  and  the  H-function  definition  of  (10)  are  compatible; 
for  practical  calculation,  however,  such  a  distinction  is  an  expensive 
luxury  and  will  not  be  maintained. 
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to  (7),  thus  yielding  the  normal  derivative  of  the  velocity  potential  in 
each  region, 


(17a)  bvCr,*)  ^  +  h K  CCts*)\ 

at 

(17b)  _  £  aj^  kjUj-r^Gss)} 

Then,  substitution  of  (17)  into  (16)  produces 

if_  -  z  +  H*r]-v  K,  r»  *>*  F*  I  « 

nee  C 

(in  which  the  prime  indicates  differentiation  with  respect  to  the  arguments- 
which  have  been  suppressed  temporarily-of  either  the  Hankei  function  or 
Legendre  polynomial).  Evaluating  this  on  the  surface  of  the  boundary  by  insert¬ 
ing  the  values  of  £.  and  £  from  (13),  and  multiplying  and  dividing  by  5*" 
leads  immediately  to 

(«>•) 


&»£>(*) /s*- 

*cj  h*o 

where,  the  finite  cylinder  wave  function  fM(s)  is  given  by 


V),  ci/eio 


(18b)  ■,‘ 

[  a,  ]a.  /*.) .  /;  r<w  »J 


which  reduces  upon  substituting  S  from  ( 12 )  to 


(19a) 

(19b) 

(19c) 


Region  I: 


Region  II 

<«> 


n 

4h  - c- 9  v&y 

Region  III: 


'ey** 


4 
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The  behavior  of  the  simplest  of  these  wave  functions,  namely,  that  of 

order  n •>  o  is  exhibited  in  Figure  2'  '  in  which  the  real  part  of  the 

wave  function  is  plotted  against  the  imaginary  part.  It 

should  be  pointed  out  that:  for  even  Y\ }  ^ „  (JT-&)  and 


X  J 

f  (*)  *  f  (if -6)  i-e.,  the  wave  functions  are  even  with  respect  to 

^  *  i  nr  i  i 

&  -  and  for  odd  (»)  =  -  ^  (IT- a)  and  (t)  =  _  f „  (ft-e) 


i.e.,  the  wave  functions  are  odd  with  respect  to  6mV/x. 

Hence,  if  the  Neumann  conditions  of  (b)  are  inserted  into  nd) 


the 


normal  derivative  of  the  velocity  potential  on  the  boundary,  it  is  readily 
seen  from  (l6)  that 


(20)  lY(fl)  = 

V  h  c 


=  ) 


3  f>(*)  I 

C  y/ihJ** 


y-  //* 


Cy/t'nJmr 


and,  therefore,  that  the  complex  coefficients  <2,,  a-h  can  be  determined, 
in  principle,  from  the  relation 


(21)  me)  -  Hx 


<**„  "A, 


(in  which  ef)^  (d)  =  )  thus  solving  the  finite  cylinder  problem, 

in  principle.  Unfortunately,  however,  it  is  not  possible  to  follow  the 
usual  procedure  of  determining  the  coefficients  of  a  general  Fourier  expansion 
and  develop  an  exact  expression  for  the  coefficient  in  terms  of  an  integral 
involving  <p*  (b)  and  tl(b)  since  different  *»<•)*  are  not .orthogonal  in 
the  0-  6  —  IT  interval^ ^ ^ .  Therefore,  a  weighted  least  squares  procedure 
which  leads  to  approximate  values  for  the  coefficients  An  and,  thus,  to  an 
approximate  solution  tor  the  finite  cylinder  problem  has  been  utilized. 

The  calculation  consists  in  approximating  the  infinite  series  Ti(B)  of 


For  this  purpose  the  to1 lowing  numerical  assignments  have  been  made 
-%a  -  r,  «=/,  a,  ;  these  assignments  are  maintained  in  all  subsequent 
numerical  calculations. 


This  is  readily  verified  after  some  numerical  analys.is  by  plotting  the 
product  of  (*)  A.  (&)  for  >"#*7  against  the  polar  angle  &  .  If  +n(e) 

is  orthogonal  to  ,  then  the  areas  under  the  real  and  the  imaginary 

parts  of  the  product  function  must  vanish  separately;  this  did  not  occur 
in  the  case  yn  =  v>  =  9,  for  example. 
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(21)  by  a  Unite  series  G(d)  of  the  form 

V 

(22a)  <o(6)  m  *2-  a.  <6  Cd) 

hmo  h  n  j 

and  requiring  that  the  weighted  average  value  of  the  square  of  the  misfit  be 
a  minimum  with  respect  to  the  choice  of  the  coefficients  Q.  on  the  boundary. 
That  is,  the  "misfit  integral" 

IT 

(22b)  t  •  J [*(•)  -  6(*)T  W(0)  M 

O 

(in  which  w(a)  is  a  suitably  chosen  weighting  function)  is  to  be  minimized 
with  respect  to  the  choice  of  by  means  of  the  requirement  that 


(22c)  „  O. 

Since,  the  are  in  general  complex  because  of  the  presence  of  the 

Hankel  functions  and  their  first  derivatives  in  the  ,  the  symmetric 

form  of  (22b)  must  be  employed  which,  when  combined  with  (22a)  and  (22c), 
yields 


(23a) 


or,  carrying  out  the  indicated  operations  and  noting  that  =  and 


-  0 


(23b) 


TT  H  IV  N 

f  s  <t>  (»)  +  Z  Z  s„*  a,*  <*„C&)  4*  (»)  l  W(»)  olfi  .  0  . 

)  I  »«*  "  Hso  m*o  ) 


By  substituting  N*(*)  -  •  Ma.  and  inverting  the  order  of  summation  and 

integration,  (23b)  becomes  the  set  of  li  y.  /  equations 

_  TT 


(2Ua) 


z  4  j»h  i"  A 
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tr 

H  (&)  W(*)  At.  even 

o  * 

which  are  the  N  relationships  required  to  calculate  the  "best"  approximate 
coefficients  at«  ... t*H  once  a  suitable  weight  function  is  selected 

The  choice  of  a  particular  weight  function  W(d)  can  not  be  made  on  an  a 
priori  analytic  basis  but  ratner  must  be  made  by  an  appeal  to  the  physics  of 

the  problem.  If  the  Riemann  integral  on  the  right-hand  side  of  (24b)  is  cast 

(7) 

temporarily  as  a  Stieltjes  integral^' ’  over  the  boundary  surface 

(25  )  /  h a  *1(6)  U t(6)  J6  ///*<*£<*)  ^  ^ 

«  s 

then,  clearly, <J\J  weights  fe)  as  some  function,  not  yet  defined,  of  the 
boundary  surface.  Consequently,  in  view  of  the  fact  that  no  preference  for 
any  portion  of  the  surface  at  the  expense  of  some  other  portion  is  inherent 
in  the  original  formulation  of  the  problem,  the  appropriate  selection  appears 
to  be  that  of  weighting  the  with  respect  to  an  element  of  surface  in 

a  manner  independent  of  the  element's  position,  that  is 

(26a)  JV  — >  d  (surface), 

or,  since  the  surface  is  a  surface  of  revolution 

(26b )  iW  -  J.TTS  ■  Sin  6  J 6  • 


or  equivalently, 

N  1 r 

(s1*)  2  »AS*ici  4>„(»)  w(6)  = 

Be®  "a  *  J 


If,  the  indicated  operations  are  performed  with  (12)  and  (15)>  then 
(27a)  dW  =  Jtt[  ± 


C6)It  should  be  noted  in  passing  that  if  the  <pJ6)  are  orthogonal  with  respect 
to  the  weighting  function  W<*>  ,  then  (2^)reduces  to  the  usual  formula  for 
the  determinating  of  the  Fourier  coefficient 

(^Admittedly  the  introduction  of  a  Stieltjes  integral  is  not  necessary  because 
of  the  presence  of  the  H-functions;  this  form  does,  however,  possess  a  certain 
aesthetic  appeal  in  view  of  the  three  sharply  defined  boundary  regions. 
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where  the  ambiguity  in  the  signs  caused  by  taking  the  root  can  be  removed 
either  by  formally  requiring  that 

(27b)  j"  JW  =  jyj1  -h  #nL 

surface 


or  by  noting  that  since  area  is  positive  and  the  cosine  and  sine  are  odd 
and  even,  respectively,  around  d«ir/4.  ,  that 


(28) 


iV  =  V/(d)  ^  =  AT 


Hence,  if  the  W(»)  of  (28)  is  substituted  into  (24b),  then  the  N/a  +1 

equations  for  the  coefficients  CL  a  ...a  become 

•  >  4  *  H 


(29a) 


where 


(29b) 


and 


(29c) 


n 

Z. “• z*: » *  >** 


r  .  |  (*,-*,)  +  -£-r  "* 

I  ♦  »  att  j  - ; - ~7*~ - 


'  uK 


TT  Z' 

^  - 


(M,  -  Ht)  *  _iL_ 

_  "  '  ^  C.V  l  A 


C.A*  * 


je 


(in  which  the  <£„(&)  of  (24b)  have  been  replaced  by  f  the  f*  (&) 
being  given  by(l8b)).  By  splitting  the  single  integrals  into  three  region 
integrals  and  making  use  of  the  fact  that  the  $»(&)  are  even  with  respect 
to  for  even  H  ,  it  is  seen  that 


(30a) 


Ax)  is) 

X  *  =  j/txJ  iZWkf*)  3»»  6  C*  »  Je 

ih  *4  * 


JS 
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and 


(30b) 


■r/x 

Cx) 


JCX) 


where  the  factor  of  AT  common  to  .!«*  and  L*  has  been  eliminated  and  where 

C»)  ti)  *  ‘  * 

the  specific  forms  of  )  and  fn(*)  are  given  by  (19a)  and  (19b),  respectively. 

Because  the  structure  of  the  integral  of  (30a)  leads  to  the 

relationship  that 


(30 


i .  =  Ci^S 

(b  4.  *•'  A 


where  cH  and  0  are  arbitrary  even  values  of  fe  and  p)  ,  only  slightly  more  than 
half  of  the  integrals  required  by  (29a)  must  be  evaluated.  Further, 
this  property  means,  since  (29a)  can  be  regarded  as  the  matrix  equation 


(32) 


<a>  <1>  -  <f>  , 


(in  which  the  complex  are  the  elements  of  the  row  vector  ,  the  complex 

^  ,  }s  are  the  elements  of  the  row  vector  and  the  complex  Cl^%)  are 

the  elements  of  the  coefficient  matrix  )  that  is  Hermitian.  Hence, 

the  inversion  of  equation  (32),  or  equivalently,  (29a)  can  be  accomplished 

/  Q  \ 

readily  by  means  of  the  technique  -  fully  described  in  Frazer'  '  -  of  triangu¬ 
lar  izing  the  matrix  followed  by  successive  backward  substitutions  of  the 

into  the  resulting  W/a,  /  equations  for  the  &,, ,  which  completes  the  approxi¬ 
mate  solution  of  the  finite  cylinder  problem. 


Numerical  Analysis  of  a  Specific  Case  -  Based  on  the  theoretical  develop¬ 
ment  of  the  preceding  section,  the  following  numerical  calculations  for  a 
three  term  series  expansion  were  performed  for  the  specific  case 

lea  -  S  Cl-*.  \  3  -  /  and  h  =  A  : 

1.  The  individual  wave  functions  f„Ci)  were  calculated  for  2, 

and  Y  ;  a  graphical  representation  of  the  calculation  for  n«o 
is  found  in  Figure  2  which  was  inserted  after  the  defining  equations 
of  (19). 

(^Frazer,  R.A.,  Duncan,  W.J.  and  Collar,  A.R.,  Elementary  Matrices,  Cambridge 
Univ.  Press,  19^6.  pp.97  et  seq. 
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2.  The  integrals  of  (30a)  and  (30b)  were  evaluated  by  calculating 
the  real  and  imaginary  parts  of  their  integrands^plotting  these 
as  a  function  of  the  polar  angle  &  and  calculating  the  areas 
under  these  curves  by  means  of  repeated  planimeter  tracings. 
These  results  are  given  in  Tables  la  and  lb: 


o 

J. 

o 

.'■7 U 

-  o/*  —  ■  A  3S~  i 

-  OtJ  +■  ASA  i. 

B 

—  •OiS‘  t 

.  at*?* 

-  ./SoS'L  eft  7 S'  i 

-.037  -037 i 

-■/30&7  -f-  oA’/’S'L 

■/f*7 

TaUe  la  -  Tine  1^. 


4 

Y*? 

B 

LA 

.4*7  -  a*/  i- 

B 

■  Oi*  j-.oaoc 

T*t>le  It,  —  TV»c  j.  *• 
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3.  The  complex  coefficients  of  the  three  term  series  expansion  of 

(29a)  were  determined  by  triangularizing  the  coefficient  matrix  <J) 
followed  by  the  successive  backward  substitutions  of  CL^  and  4^ 
into  the  resulting  three  equations  for  4,  4.^  and  CL .  The  results 
of  this  calculation  are  contained  in  Table  II: 


n 

** 

c 

•  7*  +  Ma-i 

A 

j.z  +  .7  y  i 

V- 

J./  -  .  9»A  i 

Tolo!*,  TX  —  Tt»e  cl  . 

h 


As  a  check  on  the  ability  of  these  values  to  fit  the  boundary  condition  data 
of  (6),  the  real  and  imaginary  contributions  to  the  square  root  of  the  weight¬ 
ed  "misfit"  integrand  of  (2Lb)  were  calculated  as  a  function  of  the  polar 
angle  6  ^  The  results  of  these  calculations  are  presented  graphically  in 
Figure  3  where  it  should  be  noted  that  the  extent  of  the  "misfits"  is  quite 
satisfactory  since  a  three  term  expansion  is,  admittedly,  not  sufficiently 
extensive  to  fit  the  discontinuous  data  under  consideration  with  great  "pre¬ 
cision"  . 

From  these  values  for  <Xx  and  ,  a  measure  of  the  spatial  distribu¬ 
tion  of  maximum  pressure  in  the  sound  field  is  obtained  readily  for  this 
specific  case  since,  apart  from  a  scale  factor,  ]fj^  —  h(r)l.  (10>.  These 
results  are  shown  graphically  in  Figures  4a,  4b  and  4c  for  the  three  cases: 
(i.e.,  on  the  cylinder),  p  =  and  1  -  Sd  ,  respectively. 

The  presence  of  the  pressure  "spike"  on  both  the  upper  and  lower  end  of 

the  cylinder  can  be  accounted  for  by  examination  of  Figure  5  where  the  velocity 
3 ^  (r) 

or,  equivalently,  — - -  ,  which  has  been  plotted  against  the  polar  angle, 

is  observed  to  have  a  non-zero  value  in  these  two  regions.  This  velocity 


TsT 


In  these  calculations,  the  factor  of  JLV  was  not  taken  into  account  since 
it  is  a  factor  which  could  have  been  made  to  disappear  immediately  in  (23a). 


^  ^This  f o  1  lows  imaediately  from  the  relation  PCr,  i)  m  -f  — ~rx^ 
and  the  ansatt  of  (2),  that  is,  t)  -  [iwf  *(«.-*•*  }  „  *c 

where  the  phase  ang  e  is  given  by  4  -  f^p-»  ^  ft*  y 


hence,  neglecting  the  ceaxe  factor  and  the  phased  time  dependence, 


e 


;-t 
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"spike"  follows  not  only  from  the  crudity  of  the  approximation  but  also, 
and  more  fundamentally,  from  the  fact,  which  is  emphasized  in  Figure  3> 
that  the  selection  of  the  coefficients  <*•„  was  based  not  on  minimizing  the 
unweighted  "misfit"  with  the  attendant  risk  of  having  positive  and  negative 
"misfit"  cancellation,  but  rather  was  based  on  minimizing  the  weighted  squared 
misfit  in  analogy  with  the  method  by  which  orthogonal  function  expansions 
are  matched  to  a  given  function^ 

Remarks  -  Because  of  the  extensive  numerical  analysis  required  in 
carrying  out  the  calculation  outlined  in  the  preceding  section  and  because 
a  primary  desideratum  of  the  calculation  was  its  completion  in  a  "finite" 
time,  a  less  than  ideal  number  of  terms  was  employed  in  the  series  expansion. 

If  the  "rule  of  thumb"  that  the  minimum  highest  order  Hankel  function  in  the 
series  expansion  should  be  ~  is  applied  to  this  specific  problem, 

then  /«-/£■  implies,  at  least,  a  6  term  expansion.  Such  an  effort  is 
close  to,  if  not  past,  the  economic  break-even  point  for  employment  of  an 
electronic  data  processing  system  for  this  specific  cylinder  problem  alone. 

If,  in  addition,  it  is  also  recognized  that  a  trivial  change  in  the  para¬ 
meters  4,l>  and  *4  requires  a  total  recalculation,  the  employment  of  a 
high  speed  digital  computer  becomes  mandatory. 

By  employing  an  electronic  data  processor  it  becomes  feasible  not  only 
to  perform  lengthier  calculations  but  also,  because  of  this  capability,  it 
becomes  feasible  to  postpone  to  a  later  stage  of  the  calculation  an  undesirable 
feature  which  is  imbedded  in  the  structure  of  the  current  method.  This 
limitation,  which  follows  directly  from  the  expansion  of  (l8a)  in  terms  of 

the  non-orthogonal  functions  of  (l8b),  leads  to  values  for  the  coefficients  a 

h 

which  are  dependent  on  the  number  of  terms  employed  in  the  series  expansion. 
(This  means,  simply,  that  had  a  k  term  series  expansion  been  selected  for  the 
numerical  calculation  of  the  previous  section,  say,  the  new  values  of  a.^ 
and  a.*  would  differ,  perhaps  radically,  from  those  in  Table  II  calculated 


'  Although  it  is  difficult,  if  not  impossible  to  prove,  it  seems  reasonable 
to  expect  that  the  expansion  proposed  in  (22a)  converges  in  the  mean  to  Hx 
except  for  possible  spikes  (analogous  to  the  Gibbs  phenomena)  at  the 


regional  interfaces,  as  «*•  • 
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on  the  basis  of  a  3  term  expansion.)  This  term  dependence  can  be  partially 

(12) 

offset  by  means  of  the  Gram-Schmidt  Orthonormalization  Process  which 
requires  the  use  of  a  digital  computer  for  its  effective  execution  in  this 

context . 

This  orthonormalization  process  and  the  resultant  strengthening  of  the 
finite  cylinder  calculation  is  outlined  as  follows: 

a)  From  the  relations  for  the  •f* (fl)  of  (18b)  and  (19)  and  the 
relation  for  S2-  of  (15)  construct  the  functions 


(33)  $n(fl)  =  . 

b)  From  these  non  orthogonal  functions  construct  the  ortho¬ 
normal  system  1^  for  the  finite  cylinder  (relative  to  the 

weight  function  V(*J  )  by  means  of  the  prescription 

(34a)  =  <£(a)/i!*.cs>ll 

where  ll4>0(&M|  indicates  the  weighted  norm  of  and  is 

defined  by 

r  fr 

(34b)  ||  <*,o(0)  | \  =  |  J  de  j 


and 


(35a)  >it(e)  =  Tt(t)  /iir.ooii 

in  which 


( 35b )  ?,  (e)  -  <pt  (6)  -  C <p,  u),  Hje) 


T12T 


cf.  Lindsay,  R.B.  and  Margenau,  h. 
and  Sons,  Inc.,  1936.  pp. 424-426. 


} 


Foundations  of  Physics . 


John  Wiley 
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(35c) 


(3t-a) 


(3bb) 


c) 


(37) 


d) 


(38) 


(39) 


(to) 


where  the  inner  product  is  defined  by 

TT 

( U),  *?;(«))’  I 

and 

^ra)s  V^/h^II 

in  which 

\(b)  -  ~  (m),  k /*))*,<») 

and  so  forth  until  the  required  number  of  have  been 

calculated. 

Replace  the  series  expansion  for  GC<?)  in  (22)  by 
N 

G(8)a  G(d)=  2.  Cdj 

he  O 

where  the  is  the  orthonormal  set  constructed  by  means 

of  the  Gram-Schmidt  Process  outlined  in  (3to),  (35a)  and  (36). 

Repeat  the  approximation  analysis  by  minimizing  the  "misfit" 
integral 

T  3. 

E  *  \  £/v(e)  -  e(®j]  w(eJ 

o 

and  noting  that 
TT 

J w(r)  Je  =  S- 

which  is  the  orthonormalization  condition  guaranteed  by  the 
Gram-Schmidt  Process;  then 

rlr 

which  is  term  independent. 
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e)  From  (34a),  (35b),  (36b)  and  similarly  formed  relations  for 

higher  W  ,  it  is  seen  that  the  can  be  written  as  linear 

> 

combinations  of  the  f*  in  the  form 

i 

(4la)  V[.U)  -  I.*  A;;  (»)*;<«) 

(where  it  should  be  noted  that  is  a  lower  triangular 

matrix  with  ) ;  hence,  it  follows  from  this  and 

from  (22a)  and  (37)  that 

H  Hi  /y 

(4lb)  2  k  >!;(*)  »  I  Z  =  2  *;(*) 

<■  *•  i  mo  in°  1  =  0 

or,  identifying  appropriate  terms 
H 

(42)  £t;  =  i>i  Aij 

J  t'so 

which  is  the  required  result  for  the  CLj  . 

> 

f)  Consequently,  the  term  dependence  of  the  d-i  is  introduced  in 
such  a  way  that  a  minimal  amount  of  recalculation  of  previously 
determined  CL-  S  is  needed  upon  introducing  an  increased  N  ,  say  , 

J  — — 

since  adding  A.  .  +  •••+  L*.  -  to  the  a-,  calculated  for  ri 

°  ri*  i  htij  j  H  % }  J  1 

terms  satisfies  the  requirement*  of  (42). 


Ar bitrary  Finite  Bodies  of  Revolution 


This  method  of  determining  the  velocity  potential  for  the  sound  field 
associated  with  a  finite  cylinder  is  readily  extended  to  arbitrary  finite 
bodies  of  revolution.  This  is  accomplished,  simply,  by  determining  the 
number  of  distinct  regions  and  the  extent  of  each  in  the  polar  angle  & 
for  the  cross  sectional  boundary  of  the  body,  introducing  this  number  of  H~ 
functions  such  that  each  regional  H- function  is  one  within  its  region, 
and  zero  otherwise,  using  these  regional  H- functions  to  construct  the  polar 
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equation  for  the  boundary  and,  finally,  re-doing  the  analysis  of  this  report 
by  applying  appropriate  boundary  conditions  and  introducing  different  weight 
functions  (other  than  a  uniform,  areal  weighting)  when  such  changes  are 
required  by  the  needs  of  a  specific  problem. 

NOTE  ADDED  IN  PROOF:  In  the  course  of  final  editing,  an  error  was  unearthed 
in  the  numerical  evaluation  of  the  derivative  of  certain  Legendre  Polynomials. 

A  cursory  calculation  indicates  that  the  effect  of  this  numerical  error 
should  change  the  numerical  values  cited  for  the  in  Table  II  only  moderately, 
thus  preserving  the  form  of  the  results  exhibited  in  Figures  3,  and  5. 
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Chief,  Bureau  of  Ships 
Code  320 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships 
Code  333 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships 
Code  335 

Department  of  the  Navy 
Washington  25,  D.  C. 
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Chief,  Bureau  of  Ships  1 

Code  688 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  1 

Code  689A 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  1 

Code  689B 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  1 

Code  689C 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  1 

Code  b89D 

Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  1 

Code  689E 

Department  of  the  Navy 
Washington  25,  D.  C. 

Commanding  Officer  1 

U.  S.  Naval  Ordinance  Station 
Newport,  Rhode  Island 

Officer  in  Charge  1 

Supreme  Allied  Commander  Atlantic  A.S.W.  Research  Center 
U.  S.  Naval  Base 
Norfolk  11,  Virginia 

Commanding  Officer  and  Director  2 

David  Taylor  Model  Basin 
Washington  7>  D.  C. 
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Commanding  Officer  and  Director 
U.  S.  Naval  Electronics  Laboratory 
San  Diego  Yd,  California 

Commanding  Officer 

U.  S.  Navy  Mine  Defense  Laboratory 

Panama  City,  Florida 

Commanding  Officer 

Air  Force  Cambridge  Research  Laboratory 
Air  Research  and  Development  Command 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 

Commander 

U.  S.  Naval  Air  Development  Center 
Johnsville,  Pennsylvania 

Commander 

U.  S.  Naval  Ordinance  Laboratory 
White  Oak,  Silver  Spring 
Maryland 

Commander 
Code  753 

U.  S.  Naval  Ordinance  Test  Station 

Inyoken,  China  Lake 

California 

Commander 
Code  508 

U.  S.  Naval  Ordinance  Test  Station 

Inyoken,  China  Lake 

California 

Commander 

U.  S.  Naval  Ordinance  Test  Station  Annex 
3202  East  Foothill  Boulevard 
Pasadena  8,  California 

Director 

Naval  Research  Laboratory 
Washington  25,  D.  C. 
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Director  2 

Ordinance  Research  Laboratory 
Pennsylvania  State  University 
University  Park,  Pennsylvania 

Director  2 

U.  S.  Navy  Underwater  Sound  Reference  Laboratory 
Orlando,  Florida 

Commanding  Officer  and  Director  1 

U.  S.  Naval  Training  Devise  Center 

Port  Washington 

Long  Inland,  New  York 

Superintendant  1 

U.  S'.  Navy  Post  Graduate  School 
Monterey,  California 

Committee  on  Undersea  Warfare  1 

National  Research  Council 
2101  Constitution  Avenue 
Washington,  D.  C. 

ASTIA  (TIPAA)  10 

Document  Service  Center 
Arlington  Hall  Station 
Arlington  12,  Virginia 

British  Joint  Services  Mission  1 

Bureau  of  Ships 
Code  335 

Department  of  the  Navy 
Washington  25,  D.  C. 

Canadian  Joint  Staff  1 

Bureau  of  Ships 
Code  335 

Department  of  the  Navy 
Washington  25,  D.  C. 
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